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Abstract 

The local environment and the energetic properties of one 3 He atom solved in 

bulk superfluid 4 He are studied by means of the diffusion Monte Carlo method. 

The chemical potential of the 3 He impurity is calculated with a generalized 

reweighting method which allows for a reliable estimation of this quantity. 

Results for the chemical potential, radial distribution and structure functions, 

volume-excess parameter, and effective mass are given for several pressures 

and compared with available experimental data. An overall agreement with 

experiment is obtained except for the kinetic energy of the 3 He atom which, in 

accordance with previous theoretical estimations, appears to be considerably 

larger than determinations from deep-inelastic neutron scattering. 
67.40.Yv, 02.70.Lq 
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I. INTRODUCTION 



Isotopic 3 He- 4 He mixtures have deserved theoretical and experimental interest for many 
years due to their unique properties.SH Among them one may recognize the only isotopic 
mixture which remains stable at a certain 3 He concentration down to zero temperature, 
and the only liquid system in which the two quantum statistics, bosons ( 4 He) and fermions 
( 3 He), are put together and one influences the other through the interatomic potential. As a 
result of this interplay, it has been observed both experimentally and theoretically, that the 
4 He superfluid fraction decreases and simultaneously the 4 He condensate fraction increases 
when the 3 He concentration increases. On the other hand, the 3 He momentum distribution 
in the mixture appears largely influenced by the presence of 4 He showing a considerably 
larger depletion above the Fermi momentum in comparison with pure 3 He. Experimental 
information on the 4 He condensate fraction (no) and the kinetic energy of both 4 He and 3 He 
in the mixture have been recently extracted from deep-inelastic neutron scattering.!! 4 ! These 
analysis show a large enhancement of uq with respect to pure 4 He, and 3 He kinetic energies 
very similar to the ones of pure 3 He. In contrast, all the theoretical calculations^ have 
shown only an small increment of no when the 3 He concentration (x) increases (mainly due 
to the change in the total density at a fixed pressure) and a 3 He kinetic energy appreciably 
larger. 

The maximum solubility of 3 He in 4 He is x m = 0.066 at zero pressure and presents 
a maxim value of x m = 0.095 at P ~ 10 atm. These x values are sufficiently small to 
stimulate the theoretical interest in describing microscopically the limit of zero 3 He concen- 
tration which, on the other hand, has also been experimentally analyzed and a number of 
characteristic properties are nowadays available.! From a theoretical viewpoint, this limiting 
system has been studied considering a single 3 He atom solved in bulk 4 He. The most useful 
approach in the past has been the variational method combined with the resolution of the 
hypernetted chain equations (HNC) coupledB! or not! to an Euler-Lagrange optimization 
procedure. The results obtained with this approach reproduces the energetic and struc- 
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tural properties of the system with a good accuracy but the impurity effective mass appears 
slightly underestimated.!!! The application of Monte Carlo methods, both variational0 and 
ab l nUio^ to the impurity system in order to calculate a basic property as the chemical 
potential of the impurity in the bulk (///), has been seriously hindered by the fact that /zj 
results from the difference of two energy terms of order N, N being the number of particles. 
In fact, a straightforward application of the Monte Carlo method cannot estimate /zj because 
the statistical fluctuations would mask it completely. 

In the present calculation, the reported results for /zj have been obtained using a par- 
ticular reweighting procedure suitable for the diffusion Monte Carlo (DMC) method, which 
has allowed a direct calculation of /z/ with a statistical error reduced to a manageable level. 
Using this method, we have been able to obtain reliable results for /zj that fit accurately 
the experimental data from the equilibrium up to the freezing 4 He densities. The local 
environment of the impurity, reflected in the crossed radial distribution and structure func- 
tions, has been studied by means of a pure estimatori to remove the bias associated to 
the trial wave function. We have, finally, focused our attention to the calculation of the 
impurity effective mass and its kinetic energy for several densities. As in previous quan- 
tum Monte Carlo (QMC) applications,iB the effective mass is extracted from the diffusion 
coefficient in imaginary time and, in spite of some uncertainties inherent to the extrapo- 
lated estimator used in this calculation, a reasonable agreement with recent experimental 
determinationsB'El is attained. Our results concerning the kinetic energy of the impurity, 
derived from the Hellmann-Feynman theorem, avoid the residual importance-sampling de- 
pendence and show values which are definitely larger than the experimental data, as pointed 
out previously in variational!!! and path integral Monte Carlo (PIMC)@ calculations. 

The outline of the paper is as follows. In the next section we briefly introduce the DMC 
algorithm for the impurity system and present a DMC reweighting technique that permits 
a direct estimation of arbitrarily small differences. In Sec. Ill, we present the results and 
compare them with available experimental and theoretical data. We close in Sec. IV with 
the summary and final remarks. 



II. THE DIFFUSION MONTE CARLO METHOD WITH RE WEIGHTED 

CONFIGURATIONS 

The DMC method&il allows for a very accurate description of the ground-state prop- 
erties of an interacting iV-body system. In the DMC formulation the imaginary-time 
Schrodinger equation for the function /(R, t) = \E't(R) < I ) (R, t), 

- = E "A [V, 2 /(R, t) - V,, ■ (F,(R) /(R, t))] + (E L (R) - E) f(R, t) , (1) 

is turned into a stochastic process which provides a sample of configuration points R (walk- 
ers) and weights u>(R) in a 3iV-dimensional space, whose probability distribution is given 
by /(R, t). \1/t(R) is a time-independent trial wave function that acts as an importance- 
sampling function, and $(R, t) is the exact wave function of the system. In this form, 
the Schrodinger equation appears as a diffusion-like differential equation with a diffusion, 
drift and branching terms corresponding to the first, second and third terms of the rhs 
of Eq. (H), respectively. In Eq. ([I]), E L (R) = ^ t (R) _1 ^^t(R) is the local energy, 
Fj(R) = 2\E't(R) _1 Vj\E't(R) is the quantum drift force, and Di = h 2 /(2mi) acts as the 
free- diffusion constant of the i particle. At sufficiently long imaginary times the probability 
density evolves to a stationary solution given by $o(R)^t(R) ; ^o(R) being the ground-state 
wave function from which the exact ground-state energy is obtained as the average of the 
local energy El(R). 

Let us now turn to the implementation of the reweighting method. In the DMC al- 
gorithm, the distribution probability of the walkers is modified in every single operation. 
Consider in particular the stochastic process originated by the diffusion term, which is a 
random gaussian displacement R — > R'. The new weight and distribution probability are 
w'(R') = w(R) and p'(R') = / e~ ( 4da* p(R) dR, respectively. In this stochastic process we 
can make use again of importance sampling in order to perform a modified diffusion random 
displacement. In this case, if the transition probability of going from R to R' following a 
modified diffusion process is G(R! — R), the new distribution probability is given by 
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(2) 



The statistical sample of walkers provides unchanged averaged values if one uses accordingly 
a new weight given by 



This means that a system can be studied using a variety of diffusion random laws, although 
the efficiency of the method will be related to the magnitude of the changes. In general, the 
modification has to be small enough so that the same configuration space is sampled. 

The reweighting method is specially useful in the calculation of differences between two 
almost identical systems. Performing independent samplings for both systems generates a 
global uncorrelated noise that prevents a direct measure of the difference. Assume, however, 
that given a common starting configuration point, a single drift process brings both walkers 
to new positions Ri and R 2 which are very close (in particular, separated a distance much 
smaller than the typical size of a random gaussian displacement). The configuration region 
attainable after a combined drift and diffusion process is the same, and the transition prob- 
abilities Gi(R' — Ri) and G^R' — R2) are almost equal. Equations || and |3| may be used 
then to change Gi(R' — Ri) into (^(R' — R2)- Therefore, there is no need of taking averages 
using two independent walkers for the two systems, and it may be highly preferable to use 
correlated walkers, in the sense of carrying a single random walk to obtain statistical values 
for both systems. Furthermore, notice that this technique may be applied to modify the 
diffusion process of the whole walker, i.e., all the particles of the system, or only a subset of 
it. 

The generalized reweighting method is an appropriate tool for studying the quantum 
liquid in which we are now interested. It is composed by N — 1 4 He particles and one 3 He 
atom (J) enclosed in a simulation box with periodic boundary conditions. The Hamiltonian 
of the system is 



(R'-R) 2 




(3) 
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N 




(4) 



i=i 



i<j 
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and the trial wave function ^(R) has been chosen to be of Jastrow type 



N 



* r (R)=exp(E^v)l (5) 

without distinguishing between the (4, /) and (4, 4) pairs of particles. This simplification 
in the wave function, known as average correlation approximation (ACA), has been used in 
several variational calculations!! obtaining a quite good description of the impurity prop- 
erties. In the DMC method the trial wave function acts only as a guiding wave function 
for the walkers driving them to regions where 3>o(R) is expected to be large and thus a 
particular choice, as the ACA one in the present case, does not bias the expected value 
for the ground-state energy. On the other hand, the sign problem that would emerge in a 
simulation of a finite 3 He concentration in 4 He does not appear here and an exact energy 
for the system, apart from statistical uncertainties, can be safely obtained. 

From the energetic viewpoint, the more fundamental quantity in the study of the 3 He 
impurity in 4 He is the impurity chemical potential or binding energy 

^=(H(N + I)) N+I -(H(N)) N , (6) 

both energy estimations being evaluated at fixed volume Q = N/ p. If the total number of 
particles is also conserved, and therefore one 4 He atom is substituted by the 3 He impurity, 
Hi is given by 



Hi = fi 4 + {(H((N - 1) + I)) {N _ 1)+ i - (H(N)) N ) . (7) 

We have chosen the second option in which the difference between the two energy estimations 
is much less density dependent than in Eq. (^), and moreover because it is more convenient 
if a correlated estimation of the difference is intended. The pure 4 He chemical potential /j.^ 
entering in Eq. (|7|) has been determined in a previous work using also the DMC method 
with a nice agreement with experimental dataJ^S 

The drawback of an ab initio MC estimation of fii, that has precluded such a calculation 
for years, is that an independent calculation of (H((N — 1) + -OW-ih/ an< ^ (-^(^O)iv 
lowed by its difference, produces a result completely hidden by the statistical error. In order 



to overcome this serious problem, we have directly sampled the difference by means of the 
reweighting method above introduced. Our purpose was to perform two correlated DMC 
runs, one of bulk 4 He and the other with one 3 He impurity. In this case, Eq. (|2|) has allowed 
us to use the same environment for both the 3 He atom in the impurity system and the 
equivalent 4 He atom in the pure liquid. In fact, the drift of the surrounding N — 1 particles 
in Eq. (H) is almost insensitive to the mass of the impurity, i.e., the resulting positions in 
the impurity system (R^ r _ 1 ) and in the pure phase (R^r-i) are very close. If one decides to 
change the diffusion process of the environment in the pure system with 

G^ - R-at-i) - exp (- iR ' N -l~*?- l)2 ) , (8) 

then the distribution probability of the environment, after the double process of drift and 
diffusion, is the same for the two systems. In this form, the statistical fluctuations coming 
from regions far from the impurity and its corresponding 4 He atom cancel exactly, and 
the remaining signal corresponds only to their local environment making feasible a direct 
estimation of ///. 

In addition to the impurity chemical potential ///, the knowledge of other properties as 
the crossed radial distribution function g^'^(r), the impurity effective mass and its kinetic 
energy are also relevant in a microscopic characterization of the 3 He impurity. Expectation 
values of operators O that do not commute with the Hamiltonian H are however biased 
because the probability density is \l / r(R) < I ) o(R) an d not |$ (R)| 2 - Thus, the natural ex- 
pectation values, called mixed estimators (m), have to be corrected in order to reduce or 
eliminate this systematic source of error. In the extrapolation methods^ this correction is 
approximated by 

(0) el = 2 (0) m - (0) v , (9) 

or 

{o) eq = , (io) 
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(0) v being a variational Monte Carlo (VMC) estimation. Both (0) e \ and (0) eq are accurate 
to first order in 5^, with = $o + S^f, but in general it is not enough to completely 
eliminate the influence of ^>t i n (O). In order to go beyond this approximation, we have 
used for the expectation values of coordinate operators the pure estimators following the 
methodology of Ref. O based on the future- walking strategy.il As proved in pure 4 He,0 the 



pure estimator removes all the dependence on providing results as exact as the ones for 
the total energy. 

Derivative operators as the kinetic energy cannot be evaluated with the pure estimator, 
and the extrapolation methods generate more unreliable results than in the case of 0(R). In 
a pure phase it is not a severe problem because the kinetic energy can be calculated through 
the difference E/N — V/N, V/N being the pure estimation of the potential energy. That it is 
not obviously possible in the impurity system because the total energy includes the kinetic 
energy of the medium and the one of the 3 He impurity. To overcome this difficulty and go 
to an unbiased estimation of the 3 He kinetic energy one can invoke the Hellmann-Feynman 
theorem.il It states that 

f) W 

(T,)-D,— , (11) 

E being the exact ground-state energy. We have then evaluated Tj discretizing the derivative 
dEj dDj and computing the difference in the total energy (with ADj/ Dj = 0.1 — 10% ) 
by means of the generalized reweighting method. 



III. RESULTS 

The microscopic properties of a 3 He atom immersed in bulk 4 He have been investigated 
putting it in a simulation box with N — 1 4 He atoms in such a way that the volume is 
Q = N/ p, with p the input density. In all the simulations iV = 108 particles have been used 
and the time step and population bias have been analyzed in order to remove any systematic 
error. We have also verified that for N > 100 the finite-system size introduces an error which 
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is smaller than the statistical noise, indicating that the influence of the replicas of the 3 He 
impurity implied by the use of periodic boundary conditions is negligible. The interatomic 
interaction, which does not distinguish between the two isotopes, is the HFD-B(HE) Aziz 
potential which has proved its high accuracy in a DMC calculation of the equation of 
state of superfluid 4 He at zero temperatureSS Concerning the trial wave function (|5|), the 
two-body factor proposed in Ref. ^ with the parameters optimized for pure 4 Helil has been 
considered. 

We present the results of our calculations starting with a microscopic analysis of the local 
environment of the 3 He impurity in the medium. This information is mainly contained in 
the crossed two-body radial distribution function g^'ir). In Fig. 1, mixed (short-dashed 
line) and pure (solid line) estimations of g^'ir) at densities 0.365, 0.401 and 0.424 ct~ 3 
(a = 2.556 A) are reported. In all the three densities the pure or exact results appear 
shifted to the right with respect to the mixed estimations pointing to a larger hole that is 
absolutely absent in the trial wave function. On the other hand, the height of the main peak 
in the pure <r 4 '"^(r) is slightly reduced at positive pressures and remains unchanged at the 
equilibrium density. In Fig. 2, the evolution of g^'^(r) with density is compared with the 
one shown by the pure 4 He distribution function g^'ir). Both functions show an increase 
of the localization when the density increases as well as a shift of the main peak to shorter 
interparticle distances. At a given density, the height of the main peak of <^ 4, ^(r) is smaller 
than the one of g^' 4 \r) and, what is more relevant, it appears localized to the right of the 
main peak of g^ ,4 \r) pointing manifestly to the existence of an excluded- volume region due 
to the smaller mass of the 3 He atom. The size of the excluded volume decreases when the 
density increases as one qualitatively can see comparing g^'(r) and g^'^(r) at equilibrium 
and at the highest density plotted in Fig. 2. 

Additional information on the local environment of the impurity is contained in the 
crossed static structure factor S^ 4 ' T \k), 
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N-l 



(12) 



which corresponds to the Fourier transform of g( 4 ' 7 )(r) 

S^\k) = pjdv e ikr ((? (4 '%) - 1) , (13) 

p being the density of pure 4 He. From the above definition it is easy to check that the value 
of S^\k) at the origin iM 

S^(0+) = -(l + a) (14) 

with a = v/ the quotient between the molar volume of the impurity system (v) and that 
of pure 4 He (v±). In Fig. 3, S^'^^k) is plotted in comparison with (k) — 1 at the 4 He 
equilibrium density, S^'^(k) being the pure 4 He static structure factor 

S^\k) = l + pjdr e ik - r (£ (4 ' 4) (r) - 1) . (15) 

The function S^'^(k) shown in the figure has been obtained Fourier transforming g^'^fr) 
for values k > 1 A -1 and by a direct calculation of Eq. fll^D for k < 1 A -1 . The main peak of 
S^'^(k) appears slightly depressed with respect to the one of S^'^^k) — 1 reflecting the same 
feature observed in the comparison of the radial distribution functions. Nevertheless, the 
largest differences between the two static structure functions are at low k values (k < 1 A -1 ). 
In spite of the impossibility of calculating S^' a \k) below a certain /c m i n , imposed by the use 
of a finite-size simulation box and periodic boundary conditions, if a linear extrapolation 
to k = is carried out one obtains 5^ 4 ' 4 ^(0) — 1 ~ — 1 and S^ 4 '^(0) ~ —1.3. If the latter 
is compared with Eq. (|14]) , it results a = 0.3 to be compared with the experimental value 
a expt = 0.284.0 The volume-excess parameter a decreases with pressure but this feature may 
be hardly observed in the limiting behaviour of S^'^(k) at different densities (Fig. 4). 

One of the most relevant magnitudes in the study of the impurity system is the binding 
energy of the 3 He atom in the medium or, otherwise, the chemical potential of the impurity 
pi. In Table I, we report DMC results of the pure 4 He chemical potential /14 and pi at 
three densities which correspond to the pressures also contained in the table. The results 
for the pressure and \x± reproduce the experimental data with high accuracy as pointed out 
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in Refs. p!9| , p0| . Also, in the present case, one gets a nice agreement between the calculated 
fii and the experimental dataJH the statistical uncertainties in the values of p/ being less 
than 10 %. A more exhaustive comparison between theoretical and experimental values for 
fii is displayed in Fig. 5. In the figure, two additional results are plotted: one at a pressure 
higher than 20 atm, and another located at a density smaller than the equilibrium one (po) 
which corresponds to a negative pressure of -6 atm. The solid line is a polynomial fit to the 
DMC results and has to be compared with the available experimental data of Ref. [I], a ls° 
reported in the figure. As one can see, the agreement between theory and experiment is 
excellent and a minimum in p/(p) is not observed in this region. In fact, if a minimum exists 
it is located at lower densities, even below the spinodal density of 4 He (p s = 0.264 <r -3 ).i§ 
It is worth noticing that 3 He energetically prefers to remain in the surface of liquid 4 He 
forming an Andreev state rather than penetrate in the bulk.@0 We have verifiedil that if 
the 3 He impurity is replaced by a H2 molecule there is a minimum of p/ at a density below 
Po that nearly coincides with the local density of the preferred location of H 2 in 4 He clusters 
obtained in a DMC calculation of Barnett and Whaley.il 

In ACA the chemical potential of the impurity is given byH0 

pf A = P4+(^-l)T 4 , (16) 

i. e., it can be calculated from the knowledge of properties of the pure liquid. This approx- 
imation provides upper bounds (see Table I) which, using DMC values for p 4 and T 4 , come 
close to the DMC and experimental values. 

The volume-excess parameter a may be obtained from the knowledge of p/(p), or equiv- 
alently p/(P), through the thermodynamic relation 

The values for a so obtained are reported in Table I in comparison with the experimental 
data of Ref. [TJ The agreement between a and a expt is very good at zero and intermediate 
pressures and even at high pressure where the error bar is somewhat larger. 
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Microscopic quantities which are also significant in the present study are the kinetic 
energy of the 3 He atom and the mean potential energy 3 He- 4 He (Vj). In Table II, results for 
the kinetic and potential energies for the two helium isotopes are reported at several densities. 
All of them correspond to pure estimations. In both systems, pure liquid 4 He and liquid 4 He 
with one 3 He impurity, the potential energies may be obtained using the same method that 
has been used for the radial distribution functions because they are coordinate operators. 
The pure 4 He kinetic energy simply results from the difference E / N — Vj N but that is not 
the case for Tj in the impurity system due to the coexistence of the two isotopes. Therefore, 
the kinetic energy of the impurity has been calculated using the Hellmann-Feynman theorem 
as commented in Sect. II. The ACA estimation of the partial energies of the impurity is 

T ACA = m J mi y 4 and yACA = y^ thg valueg q{ yACA bdng exp l icit l y given in Ta bl e II 

to be compared with the exact results. In going from T 4 to Tj one can see that the largest 
change is due to the difference in the mass of the two isotopes, the only effect contained in 
T ACA , and the correction due to different correlations, i. e., T ACA — T/, is in all cases less 
than 10 %. This small correction is also observed by comparing V4 and Vj. In the range of 
densities here analyzed, it is observed that Vj is always smaller than V 4 (in absolute value) 
whereas the difference T ACA - 7} is not monotonous: at P > 0, T ACA > T I but T 7 ACA < T, 
at a density 0.328 (j~ 3 (P = —6 atm). This striking behaviour can be better understood 
looking at the differences between g^cli 7 *) = 9 ( ~ A,4 K r ) an d g^'^ir) at each density. In the 
region of positive and zero pressures the main peak of g^ 4,I \r) is ever shifted to the right 
with respect to the one of g^^{r) and with a smaller localization (Fig. 2). The environment 
of the impurity may then be made equivalent to a pure 4 He liquid at a reduced density. The 
reduced density p r of the equivalent system at positive pressure can be obtained by looking 
for the density of pure 4 He at which Vj and g^ A,I \r) do correspond. If the density p r is then 
used to estimate the kinetic energy of the impurity, Tj(p r ) = 777.4/ m i ^(Pr), the results 
for Tj are the same than the ones reported in Table II. This supplies an additional test 
to our pure computation of Tj using the Hellmann-Feynman theorem. In the case of the 
equilibrium density (po = 0.365 cr~ 3 ) p r = 0.358 a -3 . At p r , we have performed an explicit 
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calculation of 4 He and have verified that g( 4 ' 4 )(r) is very much the same that </ 4,/ )(r) at 
Pq. On the other hand, at the lowest density reported in Table II (p = 0.328 cr~ 3 , P = —6 
atm) the equivalent system does not exist because the shift of the main peak of g^ 4 ' 1 '^) with 
respect to the one of g^ ,A \r) disappears and only a small derealization remains. 

There is only a previous ab initio calculation of Tj at the 4 He equilibrium density using 
PIMC and extrapolating to zero temperature.ll! Our present result for Tj, which is more 



accurate than our preliminary result of Ref. 31, is appreciably larger than the value reported 



in Ref. [12], Tj = 17.1(1) K. As a kind of closure test of our results we have calculated the mass 
dependence of Tj in order to estimate the chemical potential of the 3 He impurity through 
the relation 

(•7714 Trim) 

fij = /i 4 + / dm -!±-L . (18) 
Jmi m 

In Fig. 6, results for Tj using different masses for the impurity are displayed in comparison 
with the ACA prediction (dashed line). For simplicity, the kinetic energies Tj correspond in 
this case to mixed estimations, since at po an d for the aforementioned trial wave function 
the mixed and pure results coincide for both nij = m 4 and mj = m 3 . The PIMC result 
for m; = m 3 is also shown as an open circle. In the ACA case, if Tj in Eq. flT8] ) is 
replaced by T ACA one recovers the ACA expression for \ij (|l^) and the corresponding result 
reported in Table I, fi ACA = —2.58 K. The solid line in Fig. 6 corresponds to a fit T/(ra) = 
am + bj m, and when integrated in Eq. fllBD one obtains fij = —2.70(10) K which is 
consistent with both the experimental value and our direct estimation contained in Table I. 
As a supplementary result, it is predicted a linear departure from the ACA prediction with 
the impurity mass as is clearly manifested in Fig. 7, where the function T ACA (m) — Tj(m) 
is shown. Finally, it is worth mentioning that our results confirm and even enlarge the 
discrepancies between deep-inelastic neutron scattering determinations of the 3 He kinetic 
energy in liquid 3 He- 4 He mixturesiS (T 3 = 11 ± 3 K at P = and x = N 3 / N = 0.10) 
and all the theoretical predictions.!il@ One of the reasons that may explain this disturbing 
difference is the importance of the high-energy tails in the dynamic structure function which 
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largely influence the second energy- weighted sum rule from which the kinetic energy is 
extracted. 

We close this section with the results obtained for the impurity effective mass rrij, which 
has been recently measured with great accuracy in 3 He- 4 He mixturesBIll and also micro- 
scopically analyzed using correlated basis function (CBF) theory.0 The 3 He effective mass 
plays a relevant role in the study of 3 He- 4 He mixtures characterizing the 3 He excitations at 
low momenta. In a DMC calculation, the impurity effective mass can be obtained from the 
diffusion coefficient of the impurity in imaginary time@ 



= , im |r,(T)-r,(0)|' 
rrij r ^°° QDj r 

with Dj = h 2 / (2mj) the free-diffusion constant of the impurity. In Fig. 8, extrapolated 
estimations of rrij/ m* are reported at densities 0.365, 0.401, and 0.424 ct~ 3 . The impurity 
effective mass is extracted from a linear fit to the flat asymptotic regime of that function 
(|j~9|) which, as the figure shows, is acquired at relatively short diffusion times. The results so 
obtained are reported in Table III in comparison with the experimental determinations from 
Refs. |T4] , |r5| and the recent CBF calculation of Krotscheck et a/..0 Obviously, the experi- 
mental values are not direct measures but extrapolations to zero 3 He concentration (x) of 
determinations in 3 He- 4 He mixtures. As pointed out by Krotscheck et a/.0 a linear extrapo- 
lation, primarily used in the experimental works, is not satisfactory because the Fermi-liquid 
contributions are the most relevant in the 3 He-concentration dependence of m\ and these 
terms introduce fractional powers of x in the analytical model for rrij{x). The experimental 
values reported in Table III have been obtained using this more accurate extrapolation. 
Within the statistical errors of the DMC results, an overall agreement between our cal- 
culation and experiments is attained, with somehow a significant difference at the highest 
density due in part to the use of the extrapolated estimation (pf). On the other hand, the 



CBF results of Ref. 10 come close to the DMC and experimental results but seem to be 



slightly smaller at the densities here reported. Another CBF calculation, due to Fabrocini 
et a/.,il reported several years ago a result of m\ = 2.2 at the equilibrium density in better 
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agreement with the present DMC results. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed in this paper the most important magnitudes which characterize the 
static properties of a single 3 He atom embedded in bulk superfluid 4 He. The difficulties 
of an efficient calculation of the binding energy of the impurity in the medium, one of the 
main objectives of the present work, had prevented in the past the application of ab initio 
Monte Carlo methods to this problem. In order to overcome these difficulties, it has been 
proved that the use of reweighting techniques can be readily extended to diffusion Monte 
Carlo algorithms. This generalized reweighting method has provided reliable results for fij 
which are in excellent agreement with experimental data. 

The local environment of the 3 He atom has been explored through the calculation of 
the crossed radial distribution and static structure functions for a wide range of densities. 
The use of pure estimators for these quantities removes the uncontrolled bias, remanent in 
the approximate extrapolation methods, and shows clear evidence of an excluded volume 
region surrounding the 3 He impurity. The low k behaviour of S^'^(k) also points to the 
expected value related to the volume-excess parameter a, but a precise value for a cannot 
be estimated due to the absence of data for k < k min = 2ir / L, with L the side of the 
simulation box. Nevertheless, an independent and more precise estimation of a, through 
the pressure dependence of the chemical potential of the impurity, produces results which 
compare favorably with experimental data. 

Special attention has been devoted to an accurate estimation of the partial energies, 
potential and kinetic, of the impurity. The usual forward walking methodology does not 
apply for derivative operators, and for this reason, we have used the Hellmann-Feynman 
theorem combined with the generalized reweighting method to calculate the 3 He kinetic 
energy. The results for Tj obtained with this method show smaller differences with the 
ACA values than a previous PIMC estimate,@ with a difference T ACA — Tj which increases 
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linearly with the mass of the isotopic impurity. Our results confirm the gap between all the 
theoretical results for Tj and the much smaller 3 He kinetic energies derived from the neutron 
scattering data of Refs. |3|,£|. 

A final concern of the present work is the calculation of the 3 He effective mass through 
its diffusion coefficient in imaginary time. The results obtained show a good agreement with 
recent experimental data that slightly worsens at high pressure due probably to uncertainties 
in the MC extrapolation method used in the estimation of m*. A natural extension to the 
present work would be the calculation of the excitation energy of the 3 He impurity in liquid 
4 He, which in the limit q — > is given by h 2 q 2 / 2mj, and therefore will provide another 
method to estimate the impurity effective mass. In such a calculation, one can use DMC 
combined with the fixed-node and released-node methods, that we have already employed 
in the study of the phonon-roton spectrum in superfluid 4 He.i! Work in this direction is in 
progress. 
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TABLES 

TABLE I. Chemical potential of pure liquid 4 He (/m), chemical potential of the 3 He impurity 
(jU/), and excess- volume parameter a at several densities. The experimental data is from Ref. |]. 



p (a- 3 ) 


P(atm) p A (K) 


M ACA (K) 


W (K) 


M/ Xpt (K) 


Q ,expt 


0.365 


0. -7.27(1) 


-2.58 


-2.79(25) 


-2.785 0.284(10) 


0.284 


0.401 


10.67 -3.89(1) 


1.59 


1.38(30) 


1.42 0.200(10) 


0.199 


0.424 


20.42 -0.97(2) 


5.10 


4.73(35) 


4.83 0.176(20) 


0.165 


TABLE II. Kinetic and potential energies 


of the pure liquid 4 He and of the 3 He impurity 


immersed 


in bulk 4 He. All the ener; 


gies are in K. 






p (a' 3 ) 


T A 


v A 




Tj 


Vi 


0.328 


11.99(8) 


■19.14(6) 


15.91(8) 


17.0(6) 


-18.2(5) 


0.365 


14.32(5) 


■21.59(5) 


19.00(7) 


18.4(5) 


-21.1(5) 


0.401 


16.73(9) 


■23.88(9) 


22.20(12) 


20.5(5) 


-22.6(5) 


0.424 


18.57(8) 


■25.45(8) 


24.64(11) 


23.4(8) 


-24.7(5) 


TABLE III. 3 He impurity effective 


mass at several densities. The CBF results are from 


Ref. 0. 


p ( CT - 3 ) 


m*j 


m* cxpt (Ref. @) 


m* expt (Ref. HD 


m *CBF 


0.365 


2.20(5) 


2.18 




2.15 


2.09 


0.401 


2.36(8) 


2.44 




2.39 


2.34 


0.424 


2.72(10) 


2.64 




2.62 


2.55 
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FIGURES 

FIG. 1. Mixed (dashed line) and pure (solid line) estimations of g^ 4,I \r) at densities 0.365 
0.401 a -3 , and 0.424 cr~ 3 , from bottom to top. A vertical shift has been introduced at 0.401 
and 0.424 cr~ 3 to better visualize their differences. 
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FIG. 2. Pure liquid 4 He (dashed line) and impurity-medium (solid line) two-body radial distri- 
bution functions at densities 0.365 <r -3 , 0.401 a~ 3 , and 0.424 a" 3 , from bottom to top. A vertical 
shift has been introduced at 0.401 ct~ 3 and 0.424 a" 3 to better visualize their differences. 
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FIG. 3. Pure liquid 4 He (dashed line) and impurity-medium (solid line) static structure factor at 
the 4 He equilibrium density po = 0.365 a~ 3 . We have plotted S( 4 ' 4 )(/c) — 1 for a better comparison. 
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FIG. 5. Chemical potential of the 3 He impurity as a function of the density (full circles). The 
solid line is a polynomial fit to the DMC results. The open circles are experimental data from Ref. 

1 




0.32 0.36 0.40 0.44 



P (a -3 ) 
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FIG. 6. Kinetic energy of the impurity as a function of its mass (full circles and solid line). 
The dashed line corresponds to the ACA prediction. The open circle is the PIMC result from Ref. 
12. 




FIG. 7. Difference between the ACA prediction and the real value for the kinetic energy of the 
impurity as a function of its mass (full circles and solid line). The open circle is the PIMC result 
from Ref. 12. 




FIG. 8. The inverse of the impurity effective mass from the long-time behaviour of its diffusion 
coefficient. The solid, dashed and dotted lines correspond to densities 0.365 ct~ 3 , 0.401 cr~ 3 , and 
0.424 <r~ 3 , respectively. 
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